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ABSTRACT 


Existing  data  of  fourteen  tracer  experiments  are  analysed 
to  specify  the  pollutant  transport  and  dispersion  in  the 
north-alpine  Loisach  river  valley  for  a  variety  of  mete¬ 
orological  conditions  including  inversion  cases.  In  the 
experiments  zinc  sulfide  particles  (MMD  =  3.6  um)were  used 
as  tracers.  Particles  were  released  from  an  isolated  hill 
(300  m  abg)  at  the  valley  entrance  and  sampled  by  H-shaped 
Rotorod  devices  (located  along  the  valley  floor  and,  in  a 
few  cases,  at  some  mountain  sites  above  the  Garmisch  basin). 
Duration  of  release  was  typical  about  one  hour.  These  ex¬ 
perimental  specifics  are  briefly  outlined. 

For  unstable  lapse-rate  conditions  with  well-developed  up- 
valley  wind  (U)  the  normalized  axial  concentration  of  the 
tracer  along  the  valley  (SU/Q)  agrees  within  a  factor  of  2 
with  the  Gaussian  plume  model  (and  with  a  sigma  scheme  ap¬ 
plying  to  open  terrain  conditions) .  The  highest  centerline 
concentrations  of  SU/Q  were  observed  with  a  capping  inver¬ 
sion  layer.  In  this  case  the  measured  concentration  pattern 
was  found  in  good  agreement  with  a  Gaussian  plume  model 
assuming  perfect  reflection  from  the  inversion,  but  none 
from  the  ground. 

Tracers  were  also  detected  at  the  mountain  sites ,  in  one 
case  the  concentrations  were  considerably  higher  than  those  in 
the  Garmisch  basin.  The  transport  to  the  mountain  stations 
depends  on  the  manner  in  which  the  valley  wind  enters  the 
Garmisch  basin  and  on  the  interaction  between  the  valley 
wind  and  the  synoptic  flow  at  upper  levels. 

In  the  first  part  of  this  report  the  "straight  channel" 
model  developed  by  Gotaas  [7]  is  reviewed.  The  model  has 
been  designed  to  estimate  the  effect  of  vertical  sidewalls 
and  a  capping  inversion  layer  on  the  diffusion  from  a 
continuous  point  source.  The  model  valley  is  thought  to 
be  useful  in  case  of  reanalyzing  existing  data  or  in  plan¬ 
ning  future  experiments. 


1 .  INTRODUCTION 


Although  considerable  progress  in  understanding  the  fun¬ 
damentals  of  transport  and  dispersion  of  air  pollutants 
in  complex  terrain  has  been  achieved  in  recent  years  - 
in  particular  by  means  of  tracer  field  studies  at  various 
topographical  settings  such  as  steep-walled  canyons 
[1,  9,  20],  deep  fjord  valleys  [19]  or  broad  basins  [11]  - 
there  is  still  an  urgent  need  for  more  observational  data 
pertinent  to  terrain-related  diffusion  processes  [2] .  This 
is  most  evident  from  the  fact  that  appropriate  models  for 
impact  assessment  in  areas  of  complex  terrain  are  not  yet 
available  [19].  This  shortcoming  demonstrates  that  "ade¬ 
quate  models  for  such  conditions  will  not  be  developed 
until  the  effective  processes  have  been  better  quanti¬ 
fied"  [15]. 

The  transport  of  airborne  particles  in  a  mountain  valley 
differs  in  several  aspects  from  that  in  flat  terrain.  The 
most  striking  difference  is,  of  course,  that  the  wind- 
field  in  mountain  valleys  is  governed  to  a  high  degree 
by  thermally  driven  local  flows  which  are  often  complete¬ 
ly  decoupled  from  the  large  scale  flow.  Well  known  are 
mountain  valley  wind  circulations  and  slope  wind  regimes. 
In  those  flows  other  dispersion  mechanisms  than  small 
scale  turbulence  may  be  effective,  Fosberg  et  al.  [4], 
e.g.,  point  to  the  effects  of  organized  divergence  fields 
associated  with  the  valley  wind  and  propose  a  "divergence 
correction"  to  be  applied  to  the  conventional  Gaussian 
plume  model.  Reid  [15]  emphasizes  the  "generally  higher 
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mechanical  turbulence  levels"  leading  to  greater  dilution 
rates  near  the  source,  but  also  points  to  the  "topographic 
constraints  on  the  horizontal  dispersion  and  higher  fre¬ 
quency  of  capping  inversions"  which  may  lead  to  higher 
pollutant  concentrations  further  away. 

The  present  report  details  the  analysis  of  tracer  experi¬ 
ments  which  were  conducted  in  the  north-alpine  Loisach 
river  valley  to  investigate  the  relevant  dispersion  pro¬ 
cesses  for  a  variety  of  meteorological  conditions  inclu¬ 
ding  inversion  cases.  The  basic  data  are  summarized  in  a 
Data  Report  [17].  Since  terrain-channelling  effects  are 
important  in  mountain  valleys,  the  "straight  channel 
dispersion  model"  developed  by  Gotaas  [7]  has  been  adop¬ 
ted  to  provide  a  theoretical  frame.  This  model  has  been 
designed  "to  investigate  diffusion  in  a  valley  taking  into 
account  reflection  and  absorption  from  hillsides,  at  the 
ground  and  at  an  inversion  layer"  [7] .  Since  it  was  felt 
that  this  model  could  be  of  more  general  interest  -  e.g. 
in  case  of  reanalyzing  existing  data  or  in  planning  future 
experiments  -  the  model  is  reviewed  in  the  first  part  of 
this  report  in  detail.  A  method  is  proposed  how  the  con¬ 
centration  pattern  at  every  cross-section  of  the  valley- 
channel  can  be  obtained  without  any  cumbersome  calculation. 

2.  THE  MODEL 


2.1.  Basic  Equations 

The  model  developed  by  Gotaas  [7]  is  an  extended  Gaussian 
plume  model.  It  has  been  designed  to  estimate  the  effect 


of  side  walls  and  an  inversion  layer  on  the  diffusion  in 
a  valley  from  a  continuous  point  source.  The  reflection 
at  the  barriers  is  considered  by  optical  reflection  fac¬ 
tors  (using  the  well  known  multi  mirror  reflecting  tech¬ 
nique)  . 

Adopting  ax,  y,  z  -  coordinate  system  with  the  source  at 
P  =  (0,  0,  H) ,  the  vertical  walls  at  y  =  -  a  and  y  =  b, 
respectively,  and  the  horizontal  inversion  lid  at  z  =  h 
(Fig,  la) ,  the  diffusive  functions  S^,  and  Sz  in  the  ge¬ 
neral  relationship  for  the  concentration  S  (x,  y,  z) 

S  =  (Q/U)  Sv  St  (1) 

J  Z 


Fig.  1:  a)  Rectangular  channel  in  the  y,  z-plane 
(after  Gotaas  [7]). 

b)  The  same,  but  for  the  y*,  z*-plane. 
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emission  rate  of  the  source. 

effective  source  height  above  the  ground  surface. 

reflection  factors  of  the  ground,  the  inversion 
lid  and  the  sidewalls,  respectively. 

mean  windspeed. 
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0^  lateral  dispersion  coefficient  (standard  deviation 
of  the  crosswind  concentration  distribution) . 

oz  vertical  dispersion  coefficient  (standard  devia¬ 
tion  in  the  vertical  of  the  concentration  distri¬ 
bution)  . 

The  x-dependence  enters  functionally  since  and  oz  are 
both  increasing  functions  of  x  (see  chapter  3) . 


2.2.  Simplifications 


The  Eqs.  (1)  are  somewhat  lengthy  formulas  and  are  not 
well  suited  for  an  easy  treatment.  As  the  reflection  fac¬ 
tors  vary  between  0  and  1,  however,  a  great  deal  of  in¬ 
sight  may  already  be  gained  by  only  considering  these 
limiting  cases. 

Table  1  contains  the  simplified  formulas  for  the  horizontal 
and  vertical  diffusive  functions  for  different  combinations 
of  the  a,  0,  7-factors.  Case  I  describes  the  trapping 
between  totally  reflecting  barriers;  trapping  occurs  bet¬ 
ween  the  ground  (a  =  1)  and  the  inversion  lid  (0  =  1)  - 

see  Eqn.  (2  a)  for  the  S  -term-,  and  between  the  sidewalls 

z 

(7=1)-  see  Eqn.  (2  b)  for  the  S^-term.  Eqn.  (2  a)  is 
well  known  in  diffusion  meteorology  (e.g.  Csanady  [3])  .  In 
case  II  there  is  no  reflection  from  the  ground  (a  =  0) 
but  total  reflection  from  the  stable  layer  (0  =  1) ,  see 
Eqn.  (3  a).  Case  III  is  the  widely  used  Gaussian  plume 
model  (e.g.  Stern  et  al.  [21]),  with  total  reflection 
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Table  1 ;  Simplified  formulas  for  the  diffusive  functions 
S  and  S  for  different  combinations  of  the  re- 

y  2 

flection  factors  a,  y-  For  details  see  text. 
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from  the  ground  (a  =  1)  and  no  reflection  from  the  stable 
layer  (equivalent  to  case  of  absent  diffusion  lid) ,  see 
Eqn.  (4a).  In  case  IV  all  reflection  factors  are  zero, 
and  thus,  no  reflection  takes  place,  see  Eqs.  (5a)  and 
(4b) ,  respectively. 

2.2.1.  Some  Often  Used  Formulas 


For  many  practical  purposes  the  ground  level  concentration 
is  of  special  interest.  A  convenient  measure  is  the  wind- 
speed-normalized  relative  concentration  SU/Q  which  is  sim¬ 
ply  given  by  SU/Q  =  S^Sz  (Eqn.  1). 

When  the  Eqs.  (3a,  4a,  5a)  for  the  S  -term  are  combined 

z 

with  the  Eqn.  (4b)  for  the  S^-term  (Table  1)  one  gets  in 
case  of  the  axial  (y  =  0)  concentration  for  the  ground 
level  (z  =  0)  the  following  relations: 

a)  Total  reflection  at  the  inversion  lid  (3  =  1),  no  re¬ 
flection  at  the  ground  (a  =  0) : 

-1  r  _  0.  _  (Hi.  1 

SU/Q  -  (zir 6j  [e  ^  +  c  J  (6, 

b)  Total  reflection  at  the  ground  (a  =  1),  absence  of  an 
inversion  lid  (Gaussian  plume  model) : 

SU/Q  *  (tc  g3  gz)  e 


i' 


(7a) 
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c)  No  reflection  at  the  ■-'round  (a  =  0)  ,  absence  of  an  in¬ 
version  lid  (Gaussian  absorption  formula) : 

SU/Q  ‘(pc  &,  Gz) 1  e  (7bi 

2.3.  Dimensionless  Equations 

The  Eqs.  (2a)  and  (2b)  in  Table  1  describe  the  vertical 
and  horizontal  variation  of  the  plume  concentration  in 
a  straight  channel  with  totally  reflecting  barriers 
(a  =  0  =  y  =  1 ) ;  the  dependence  on  the  x-direction  en¬ 
ters  through  the  semi-empirical  relationships  o^(x)  and 
a  (;•)  ,  see  chapter  3.  A  few  steps  are  needed  to  bring 
this  set  in  a  more  tractable  form.  In  order  to  derive  a 
general  non-dimensional  relationship  the  vertical  and 
horizontal  diffusive  functions  are  treated  separately. 

2.3.1.  Vertical  Diffusive  Function 

Using  the  height  of  the  inversion  layer  h  as  scaling 
height  the  following  dimensionless  quantities  are  ob¬ 
tained: 

z*  =  z/h,  H*  =  H/h,  a  *  =  a  / h.  (8) 

z  z 

Note  that  both  the  height  coordinate  z*  and  the  source 
position  parameter  H*  vary  between  0  and  1  (see  also 
Fig.  1b) .  With  this  the  dimensionless  vertical  diffusive 
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function  S  *  is  related  to 
z 

Sz 

S  *  =  h  S  ,  (5) 

4  z 


where 


n*  *  qo 


(z*-  h\Znf 


t  C 


(10) 


The  case  of  no  reflection  at  the  inversion  lid  (0  =  0) 
leads  to  the  Gasussian  plume  model: 


2.3.2.  Horizontal  Diffusive  Function 

In  an  analogous  way  the  width  of  the  valley  B  (=  a  +  b) 
is  used  to  get  the  other  dimensionless  quantities: 

y*  =  (y  +  a)/B,  a*  =  a/B,  oy*  =  o  / B  (12) 

It  follows,  that  y*  =  0  for  y  =  -  a  and  y*  =  1  for 
y  =  b,  and,  therefore,  the  lateral  dimensionless  coor¬ 
dinate  y*  varies  in  the  range  0  £  y*  s  1 (see  also  Fig. 1b). 
The  source  position  parameter  a*  varies  also  between  0 


and  1 ,  in  case  of  the  source  position  at  the  center  of 
the  valley  (a  =  b  =  B/2)  it  follows  a*  =  0.5. 

The  dimensionless  horizontal  diffusive  function  S  *  is 

Y 

related  to  S  by 

y 


S  *  =  B  s  , 
y  y 


where 


00  r  ('4*~Q?+Zr>) 

+e‘ 


s3  “  isrqi.  le 


The  case  of  no  reflection  at  the  sidewalls  (y 
to  the  Gaussian  plume  model: 

c*  a  — 1 _  c“  2*^ 

j  T'SkgJ 


=  0)  leads 


2.3.3.  General  Relationship 


With  the  aid  of  the  Eqs.  (10)  and  (14)  for  S  *  and  S  * 

z  y 

the  concentration  S  can  be  expressed  as 


S/ST ,  =  S  *  .  S  * 

'  Lim  y  z 


wher^ 


0 


12 


ST .  is  the  far  field  concentration  in  the  channel,  the 
Lim 

concentration  at  great  distance  will  not  fall  below  this 
limiting  value. 


Since  Sy*  and  Sz*  show  the  same  functional  dependence, 
i.e.,  Sy*  *  F*  (o  *,  y*,  a*)  and  Sz*  =  F*  (oz*,  z*,  H*) , 
a  general  non-dimensional  diffusive  function  F*  may  be 
expressed  as: 


,  00  (£*-$%  2n) 


(e* +  f+2n) 
2  i 


ru-«* 


In  case  of  F*  *  Sy*  the  general  parameters  are  identified 

as  a*  =  a  *,  1*  -  y*,  C*  =  a*,  and  in  case  of  F*  =  S  * 

y  2 

as  a*  *  a  *,  1*  *  z*,  5*  *  H*.  Thus,  the  concentration 
z 

ratio  S/SLim  can  be  written: 


S/SLim  *  F * ( oy * ,  y*,  a*)  •  F*(oz*,  z*,  H*) 


In  order  to  enable  a  quick  estimate  of  s/sLim  according 
to  Eqn.  (19),  the  functional  dependence  of  F*  on  the  ge¬ 


neral  parameters  a*,  1*  and  5*  (Eqn.  18)  has  been  tabu¬ 
lated  (see  Table  B1  in  appendix  B) .  The  repeated  appli¬ 
cation  of  the  tabulated  set  -  both  for  the  parameters 
az*,  z*,  H*  and  for  oy *,  y* ,  a*  -  gives  s/sLim  at  anY 
desired  location  of  the  channel  in  a  very  easy  way  (an 
example  is  given  in  chapter  4) .  Of  course,  this  method 
presupposes  that  the  functional  dependence  of  the  disper¬ 
sion  coefficients  (a.,  o  or  o,  *,  0  *)  on  stability  and 

y  z  y  z 

distance  is  known.  This  problem  is  postponed  to  chapter  3, 
the  next  step  is  a  short  discussion  on  the  characteristics 
of  the  diffusive  functions. 


0 


2.3.4.  Characteristics  of  the  Diffusive  Functions 


The  functional  behaviour  of  F*  (a*,  1*,  C*)  for  a  variety 

of  specified  conditions  is  to  be  seen  from  the  Figures 

2-4.  Note  that  F*  =  S^*  in  case  of  o*  =  CTy*»  1*  =  Y*» 

5*  =  a*  or  F*  =  S  *  for  o*  -  a  * ,  1*  -  z* ,  c*  =  H* , 

z  z 

respectively. 

Fig,  2  shows  the  variation  of  F*  with  distance  a*  at 

various  locations  1*  (parameter) ,  the  source  position  is 

fixed  at  *  0.5.  Hence,  in  case  of  F*  «  Sy*  the  diagram 

depicts  the  cross-section  distribution  (0  £  y*  £  1)  at 

the  distance  for  a  source  located  at  the  centre  of 

the  valley  (a*  =  0.5).  For  F*  =  Sz*  it  shows  the  vertical 

distribution  (0  £  z*  £  1)  at  the  distance  a  *  for  a  source 

z 

located  midway  between  the  ground  and  the  inversion  lid 
(H*  *  0.5).  In  either  case  it  is  seen  that  F*  *  1  for 
o*  5  0.5. 

Fig.  3  shows  the  variation  of  F*  with  distance  a*  at  the 
location  1*  =  0  for  different  source  positions  C*  (para¬ 
meter)  .  Hence,  this  diagram  depicts  the  variations  along 
the  wall  (y*  *  0)  and  at  the  ground  (z*  *  0) ,  respective¬ 
ly.  It  is  seen  that  F*  =  1  for  a*  2  1 . 

Fig.  4  reflects  the  case  1*  =  c*  for  different  source 
positions  (parameter) .  In  case  of  F*  =  S  *  the  variation 
refers  to  the  source  height  (z*  =  H*) ;  the  curve  z*  -  0, 
e.g.,  describes  the  variation  of  S  *  at  ground  level  for 

Z 

a  ground  level  source.  It  is  seen  that  F*  approaches  the 
limit  1  in  the  range  0.5  <  o*  £  1. 
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Fig.  4:  Variation  of  diffusive  function  F*  with  non- 

dimensional  distance  o*  at  the  location  1*  =  t* 
according  to  Eqn.  (18)  .  Parameter  is  the  source 
position  C* . 

A  comparison  was  also  made  between  the  multi  reflection 

formula  Sz*  (Eqn.  10)  and  the  respective  formula  for  the 

Gaussian  plume  model  S  *  (Eqn.  11)  for  the  ground  level 

2  ~ 
z*  ~  0.  Fig.  5  shows  the  result,  i.e.  the  ratio  S  */S  * 

Z  Z 

as  a  function  of  c?z*  for  different  source  positions  H* 
(parameter) .  Since  S  *  =  1  for  a  *  £  1  (Fig.  3) ,  the 

Z  Z 

curves  for  a z*  £  1  simply  reflect  the  relationship 

1/S,*  *  V'n/2'  a  *  exp  (H*2/2a  *2)  .  For  a  *  <  1  the  ratio 
z  z  z  z 

varies  between  1  and  2.  The  factor  2  is  found  for  H*  ■  1 , 
i.e.  with  the  source  at  the  height  of  the  upper  lid,  in 
this  case  the  ground  level  concentration  is  doubled  every 

A 

where.  Similar  diagrams  may  be  obtained  for  S  */S  *. 


5:  The  ratio  S  */S  *  at  ground  level  (z*  =  0)  as  a 
—  z  z 

function  of  a  *.  Parameter  is  the  source  position 


3.  FORMULAS  FOR  THE  DISPERSION  COEFFICIENTS 


In  order  to  treat  specific  problems,  the  dependence  of 

a  ,  a  (or  o  *,  a  *)  on  stability  and  distance  must  be 
y  z  y  z 

known.  Since  there  are  several  sigma  schemes  in  use  (see 
[5,  6,  19]),  the  reader  may  apply  a  scheme  of  his  own 
choice.  In  the  following  we  rely  on  the  scheme  proposed 
by  Nowicki  (13],  which  is  based  on  a  comprehensive  ana¬ 
lysis  of  numerous  diffusion  measurements.  The  "universal" 
relationships  for  o  and  oz  not  only  consider  the  depen¬ 
dence  on  atmospheric  stability  and  distance  x  but  also 
on  height  and  aerodynamic  roughness  zQ.  The  classification 
of  stability  is  done  via  the  exponent  m  in  the  power  law 
of  the  vertical  wind  profile.  For  the  mean  values  between 
zq  and  H,  i.e.  for  the  coefficients  to  be  used  directly 
in  the  diffusion  formulas,  Nowicki  suggests  the  following 
relationships : 


oy  =  0 . 08 ( 6ra~° ' 3  +  1  -  lnH/zo)X0'367 (2'5_m) 
o  =  0.38m1 *3(8. 7  -  InH/z  )  X1  * 55  exp  (_2'35m) 


(20) 


(21) 


With  a  "roughness  function"  R  defined  by 


R  =  0.08(8.7  -  InH/z  ) , 


(22) 


the  power  laws  for  oy  and  oz  may  be  expressed  as 


o  =  Fx 

y 


(23) 


with 


0 


20 


Table  2  shows  how  the  stability  parameter  m  and  the  in¬ 
dependent  functions  Fq,  f,  Gq,  g  vary  within  the  widely 
used  diffusion  categorization  scheme  discussed  by  Pas- 
quill  [14]  and  Turner  [22]. 

With  this  frame  the  a-power  laws  (Eqs.  23,  24)  are  rea¬ 
dily  obtained:  for  a  special  problem  -  with  known  H  and 
zq  -  it  is  only  necessary  to  calculate  the  roughness 
function  R  (Eqn.  22)  to  get  F  (Eqn.  23a)  and  G  (Eqn.  24a) 
in  a  very  easy  way;  the  exponents  f  and  g  in  the  power 
laws  may  directly  be  taken  from  Table  2  for  the  respec¬ 
tive  stability  class  (an  example  is  given  in  the  next 
chapter) . 

4.  EXAMPLE  OF  DIFFUSION  SPREAD  IN  THE  STRAIGHT  CHANNEL 
MODEL 


With  the  availability  of  the  o  <x)  and  a  (x)  power  laws 

y  z 

an  estimate  of  the  concentration  pattern  s/sLim  (Eqn.  19) 
may  be  obtained  for  any  given  channel  configuration.  For 
practical  purposes  it  often  is  sufficient  to  have  some 
knowledge  of  the  concentration  distribution  at  a  few 
selected  cross-sections.  This  can  be  accomplished  very 
rapidly  and  without  any  cumbersome  calculations.  How  to 
proceed  in  this  case  will  be  outlined  in  this  section 
with  an  example. 

Starting  point  is  the  ground  level  concentration  pattern 
for  a  channel  of  width  B  =  2000  m  and  height  h  =  500  m, 
the  point  source  is  at  H  =  100  m  and  its  distance  from 


y,  m 
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Fig.  6;  Ground  level  concentrations  S(x,  y ,  0)  within  the 
’straight  channel'  downwind  of  a  continuous  point  source 
at  height  H  =  100  m  and  with  a  =  400  m  (lateral  distance 
from  the  wall) .  Numbers  at  the  grid  points  denote  S  in  10“ 
particles/m3.  For  model  input  parameters  see  Table  3. 
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the  wall  is  a  =  400  m  (see  Fig.  6) .  Numbers  at  the  grid 
points  (mesh  widths  in  lateral  (y)  and  longitudinal  (x) 

directions  are  100  m  and  200  m,  respectively)  denote  the 

3  3  4 

concentration  S  in  10  particles/m  .  For  S  £  10  par¬ 
ticles/m3  isolines  are  drawn  also.  With  an  emission  rate 
Q  =  1.3  *  101  particles/s,  a  mean  windspeed  of  U  =  5  m/s 

and  the  channel  dimensions,  the  far  field  concentration 

3  3 

£>Lim  (Eqn.  17)  amounts  to  2.6  *  10  particles/m  .  Fig.  6 
shows  how  this  limiting  value  is  approached  from  both 
channel  sides  with  growing  distance  (change  from  inhomo¬ 
geneous  to  homogeneous  conditions) . 


A  scheme  has  been  prepared  to  allow  the  estimate  of  the 
concentration  at  any  desired  distance  x  (Table  3) .  The 
model  input  includes  the  specification  of  the  diffusion 
meteorology,  the  channel  configuration  and  the  source 
parameters;  no  further  comment  is  necessary  on  the  deter¬ 
mination  of  SLim,  the  dimensionless  parameters  [z* ,  H*, 
a*)  and  the  specification  of  the  factors  in  the  power 
laws  for  Oy,  o z«  For  illustration  the  parameters  of  the 
chosen  example  (Fig.  6)  are  listed  in  the  scheme. 


The  lower  half  of  the  table  shows  the  cross-section 

distribution  of  S/ST  .  and  S  at  the  distance  x  =  3000  m 

Lim 


(example) .  Since  S/S 


C  * 

bz  1 


it  is  only  necessary 


=  c  * 

' Lim  by 

to  use  Table  B1 ,  to  pick  out  the  respective  values  for 
S  *  and  S  *,  and  to  do  some  little  desk  calculator  multi- 

y  z 

plication  (in  case  of  S^,*  it  holds:  F*  =  Sy*'  ~  a*  =  0.2, 
a*  =  a  *  ‘ 


of  S 

y 

range 


=0.16  (Table  B1).  The  cross-section  distribution 
*  is  tabulated  for  1*  =  y*  ("upper  scale")  in  the 
0  i  y*  i  1 ,  these  values  are  listed  in  Table  3) . 
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Table  3:  Scheme  of  input  parameters  for  the  straight  channel 
model  with  an  example.  For  details  see  text. 


Stability  class 

A-F 

c 

Mean  wind  speed 

U 

5  m/s 

Height  of  the  stable  layer 

h 

500  m 

width  of  the  valley 

B 

2000  m 

Aerodynamic  roughness 

0.20  m 

SOURCE  PARAMETERS 

a)  Strength 

0 

l^lQ10  particles/s 

t»  Height 

H 

100  m 

c)  Lateral  distance  from  the  wall 

a 

400  m 

Farfleld  concentration  S^  ■  fl/(U-B‘h) 

sLlm 

2.6*10^  partldes/m3 

Immlsslon  level  (z  -  0:  ground) 

z 

0  m 

Z*»  z/h 

z* 

0 

H*-  H/h 

H* 

0.20 

a*-  a/B 

a* 

0.20 

Specification  of  the  power  laws  for  the 

dispersion  coefficients  toy  -  Fxf,  o2  -  Sx3> 

H/Zo 

H/Zq 

500 

R  -  0.08  (8.7-lnH/Zo) 

R 

0.1388 

F0  (Table  2) 

Fo 

0.1666 

So  (Table  2) 

So 

0.5710 

F  -  F0>  R 

F 

0.565 

G  -  6q  •  R 

6 

0.114 

f  (Table  2) 

f 

0.846 

g  (Table  2) 

g 

0.978 

DISTANCE 

X  (a) 

V 

3Z* 

V 

5000 

0.1S 

0.57 

1.32 

eras  -  SECTION  0*y'*l 


0.0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

V 

2.28 

2.48 

2.50 

2.07 

1.14 

0.45 

0.11 

0.02 

0 

0 

0 

*\im 

5.0 

3.3 

3.4 

2.7 

1.5 

0.6 

0.2 

0 

0 

0 

0 

( 10 3  Bertie!**/*3 

7.8 

8.3 

8.9 

7.1 

3.9 

1.5 

0.4 

0.1 

0 

0 

0 

* 
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Any  other  configuration  may  be  treated  in  an  analogous 
manner.  This  procedure  may  be  very  useful  in  case  of  a 
reanalysis  of  existing  data  or  in  case  of  planning  strate¬ 
gies  for  future  experiments. 

5.  DIFFUSION  IN  A  REAL  VALLEY 

5.1.  General  Remarks 

The  fluorescent  particle  (FP)  tracer  experiments  in  the 
Loisach  river  valley  area  were  conducted  to  gain  basic 
insights  into  the  mountain  diffusion  meteorology  and  to 
test  the  applicability  of  the  conventional  model  frame¬ 
work  to  the  complex  conditions  in  a  deep  mountain  valley. 
The  details  of  the  fourteen  experiments,  which  took  place 
from  May  1975  to  July  1976  under  a  variety  of  meteorolo¬ 
gical  conditions  including  inversion  cases,  are  summari¬ 
zed  in  a  Data  Report  (prepared  for  the  European  Research 
Office,  United  States  Army,  London,  see  Reiter  et  al. 

[17])  . 

Since  the  subsequent  analysis  relies  on  these  data  some 
background  information  on  the  topography  and  the  experi¬ 
mental  procedures  will  be  repeated  here  for  the  reader's 
convenience . 

5.2.  Site  Description  and  Measuring  Facilities 

The  Loisach  river  valley,  which  is  well-known  for  its 
distinct  diurnal  wind  system  [16,  18]  -  with  daytime 
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north-eastern  (NE)  up-valley  winds  and  nighttime  south¬ 
western  (SW)  down-valley  winds  -,  is  about  20  km  long, 

2  km  wide  and,  in  its  northern  part,  1  km  deep  (Fig.  7) . 

The  northern  part  of  the  valley  is  SSW  -  NNE  oriented  and 
enters  the  'Murnauer  Moos*  fen  or  the  Bavarian  pre-alpine 
region  in  a  f unr.el-shaped  way. 

The  other  section  of  the  valley,  the  Garmisch  basin,  runs 
from  WSW  to  ENE  and  thus  shows  a  considerably  deviating 
direction.  In  the  south  it  is  surmounted  by  several  ranges 
of  the  Wetterstein  massif  with  the  Zugspitze  (t3000  m  a.s.l. 
being  its  highest  peak.  Since  the  main  ridge  raises  to 
2600  m  height  or  almost  2000  m  above  the  valley  floor, 
the  southern  ranges  are  by  far  the  highest  of  all  surroun¬ 
ding  mountain  chains  including  those  of  the  Kramer  complex 
in  the  northwest. 

The  walls  of  the  main  valley  are  forested  up  to  the  timber- 
line  at  about  1700  m  a.s.l.;  the  sloping,  however,  varies 
considerably  from  place  to  place,  only  the  eastern  flank 
(Estergebirge)  of  the  northern  part  shows  a  fairly  homo¬ 
geneous  structure  with  an  inclination  of  approximately 
30°  to  a  height  of  1300  m  above  the  river. 

The  nature  of  the  valley  floor  is  characterized  by  mea¬ 
dows,  small  forests  and  urban  districts  marking  this  area 
as  one  of  considerable  inhomogeneous  aerodynamic  rough¬ 
ness. 

The  plan  to  accomplish  diffusion  measurements  in  this 
area  has  been  considerably  promoted  by  the  existence 


Fig.  7:  Map  of  the  Loisach  River  area  with  contour-lines 
(m)  drawn  in  100  m  intervals.  Tracer  was  released 
at  the  Hdhenberg  mountain  'H'.  Samplers  were  lo¬ 
cated  at  the  valley  floor  and  at  mountain  sites 
(Wank  peak  and  sites  numbered  I  -  VI) .  The  Insti¬ 
tute  is  indicated  by  the  letter  *1'.  Dashed  lines: 
Rivers  or  creeks.  Dotted:  ridge  lines. 
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of  an  isolated  hill  (300  m  abg)  in  the  immediate  vicinity 
of  the  valley  entrance.  This  hill,  being  a  unique  topo¬ 
graphical  feature,  invites  to  release  the  tracer  from 
its  top  (the  'release'  mountain  is  indicated  by  an  'H' 
in  Fig.  7) .  With  the  tracer  released  at  the  valley  en¬ 
trance  our  primary  objective  has  been  to  investigate  the 
aerosol  transport  along  and  across  the  valley  under  a 
variety  of  characteristic,  but  different  meteorological 
up-valley  wind  conditions. 

Generally,  most  samplers  were  installed  at  various  down¬ 
wind  locations  at  the  valley  floor,  in  several  cases, 
however,  some  few  devices  were  also  run  at  selected  moun¬ 
tain  sites  (Wank  peak  and  sites  labeled  by  roman  numbers 
(I  -  VI)  in  Fig.  7)  . 

For  each  experiment  comprehensive  meteorological  informa¬ 
tion  was  provided:  i)  by  the  permanent  meteorological 
measuring  facilities  at  the  Institute  (indicated  by  an 
"I"  in  Fig.  7)  and  the  surrounding  high  mountain  observa¬ 
tories  Wank  and  Zugspitze;  ii)  by  special  pibal  tracking 
(windfield)  and  radiosonde  ascents  (temperature)  at  seve¬ 
ral  locations  in  the  valley  prior  to,  during  and  after 
each  experiment.  Cloud  cover,  radiation  conditions  and 
other  relevant  parameters  were  also  included  to  gain 
further  insight  into  the  diffusion  meteorology. 

5.3.  Experimental  Techniques 


The  tracers  were  zinc  sulfide  fluorescent  particles  (FP) 
from  the  United  States  Radium  Corporation  (USCR) .  The 
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mass  median  diameter  MMD  was  3. 64m  and  the  particle  den¬ 
sity  PPG  (particles  per  gram)  was  0.91  *  lO1^. 

The  dissemination  of  the  aerosol  was  accomplished  by  a 
Metronics  Model  8  Blower  Generator  of  the  series  widely 
used  in  the  field  (Leighton  et  al.  [12]).  For  all  expe¬ 
riments  an  emission  rate  Q  -  1.3*10**^  particles/s  was 
used,  the  duration  of  emission  varied  between  40  and  60 
min.  This  mode  of  acting  assured  sufficient  coverage  in 
all  cases  to  apply  the  FP-technique  successfully. 

Tracer  samples  were  collected  by  H-shaped  Rotorod  samplers. 
These  were,  however,  no  Metronics  fabricated  devices,  but 
in  fact  the  Metronics  standard  type  (as  described  by  Gri- 
nell  et  al.  [8],  or  Leighton  et  al.  [12])  was  reproduced 
by  our  laboratory,  with  a  total  of  20  devices.  Considering 
the  Rotorod  efficiency  of  about  65%  (for  the  particle  size 
range  used  in  the  experiments  and  for  rods  coated  accor¬ 
ding  to  the  standard  procedure  [8]),  the  true  sampling 
rate  turned  out  to  be  26.9  1/min. 

Before  each  experiment  the  collector  arms  of  the  Rotorods 
were  "manually  coated"  with  special  silicone  grease  ac¬ 
cording  to  the  recommended  standard  procedure  [8] . 

During  the  experiment  all  samplers  were  fixed  to  metal 
posts  at  approximately  one  meter  above  the  ground,  as  is 
common  practice  in  comparable  field  trials  (e.g.  Archuleta 
et  al .  [ 1 ] ) . 
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The  samplers  were  operated  on  specially  designed  9-volt 
d.c.  battery  systems  providing  constant  rotation  speeds 
(with  a  constancy  better  than  that  of  the  standard  ver¬ 
sion  (±2%)  during  a  several  hours  run)  . 

The  samplers  were  energized  just  prior  to  a  release.  After 
cloud  passage  the  period  of  operation  was  'held  to  a  mini¬ 
mum  in  order  to  avoid  obscuration  of  FP  by  atmospheric 
particulates  deposited  after  cloud  passage'  (as  has  been 
recommended  by  Leighton  et  al.  [12]). 

After  each  experiment  the  particles  on  the  collector  rods 
were  counted  by  means  of  a  Zeiss  microscope  of  magnifica¬ 
tion  160  x  (10  x  eyepiece  and  16  x  objective  of  0.35  N . A. ) 
with  incident  UV  light  (to  excite  the  fluorescence) . 

With  the  particle  counts  denoted  by  DT  -  where  the  sam¬ 
pling  time  t  (min)  is  indicated  by  the  index  -  the  mean 
particle  concentration  Sx  was  determined  by 


with  Fr  =  26.9  lmin  1  being  the  true  sampling  rate. 

The  concentration  was  converted  to  a  S^q  10  min  average 
by  the  power  law  recommended  by  Turner  [23]: 

S1Q  =  ST(x/10)0*2.  (26) 

In  case  of  t  =  60  min,  the  S,n  values  were  multiplied  by 

b  U 

1.43,  a  conversion  factor  well  known  in  diffusion  meteoro¬ 
logy. 
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5.4.  Survey  on  the  Experimental  Specifics 


A  survey  on  the  experimental  specifics  -  release  data, 
meteorological  conditions,  number  of  samplers  at  diffe 
rent  areas  of  interest  -  is  given  in  Table  4. 

Table  4 :  Survey  on  the  experimental  specifics  of  the 
Loisach  river  valley  diffusion  experiments. 


EXPERIMENT 

Meteorological 

conditions 

Njnber  of  samplers 
at 

Nunber 

Date 

Time 

(CET) 

Duration 

(min) 

Stability 

class 

Wind 

sceed 

(ms*1) 

Northern 
oart  of 
tne  valley 

Gamuscl) 

basin 

Mountain 

sites 

1 

13  May  75 

12:35 

60 

D 

3.0 

18 

1 

l 

2 

26  Jun  75 

11:00 

60 

C(B) 

6.0 

18 

1 

l 

3 

7  Jul  75 

11:10 

60 

B 

5.5 

18 

1 

l 

4 

9  Jul  75 

11:30 

60 

CCD) 

4.5 

13 

5 

- 

5 

23  Jul  75 

12:04 

60 

B 

6.0 

19 

1 

l 

6 

28  Jul  75 

12:00 

40 

C(B) 

6.5 

wm 

6  Aug  75 

11:30 

40 

CCD) 

6.0 

1 

l 

■a 

13  Aug  75 

12:00 

40 

C 

5.0 

■a 

9 

U  Nov  75 

12:45 

40 

5.5 

20 

- 

- 

10 

16  Dez  75 

13:00 

40 

J 1  I 

20 

- 

- 

U 

8  Mar  76 

11:30 

60 

5.0 

20 

- 

- 

12 

14  ADr  76 

10:15 

45 

- 

20 

- 

- 

13 

28  Jun  76 

11:00 

45 

8(0 

6.0 

6 

9 

5 

14 

7  Jul  76 

10:30 

60 

B 

7.0 

8 

9 

3 
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The  stability  class  was  determined  by  the  most  widely 
used  diffusion  categorization  scheme  discussed  by 
Pasquill  [14]  and  Turner  [22],  and  the  mean  flow  was  spe¬ 
cified  by  an  average  wind  speed  between  ground  level  and 
300  m  height  (source  level)  deduced  from  the  pibal  measu¬ 
rements.  According  to  this,  the  stability  ranged  between 
B  and  D  categories,  and  the  wind  speed  varied  between  3 
and  7  m/s.  Most  (10)  experiments  were  conducted  during 
the  summertime  with  we 11 -developed  up-valley  winds,  where¬ 
as  the  remaining  four  experiments  represent  winter/spring 
cases  with  partly  complex  meteorological  conditions  (in¬ 
version  structures  and  in  one  case  (No.  12)  unsteady 
winds)  . 

5.5.  Analysis 


The  results  of  the  FP  tracer  sampling  program  are  discus¬ 
sed  in  terms  of  the  wind-speed-normalized  relative  con¬ 
centration  SU/Q  which  appears  to  be  the  most  convenient 
entry  when  comparing  dilution  rates  of  different  experi¬ 
ments.  A  comparison  with  theory  was  also  made  whenever 
it  seemed  appropriate.  Special  aspects,  e.g.  the  release 
of  the  tracer  material  above  a  persistent,  lifted  ground 
based  inversion,  are  treated  in  a  separate  section  (5.5.3.). 

5.5.1.  Inversion  Lid  Conditions 


A  fairly  good  example  for  the  diffusion  during  an  inver¬ 
sion  with  the  source  below  the  inversion  base  is  the  Ex¬ 
periment  No.  12  (8  March  1976).  In  this  case  the  highest 
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absolute  concentrations  of  all  experiments  were  observed, 
and  also  from  this  point  of  view  a  more  detailed  discus¬ 
sion  seems  to  be  appropriate. 

Fig.  8  shows  the  particle  counts  at  the  valley  floor  to¬ 
gether  with  the  plume  centerline.  The  distribution  appears 
remarkably  well  organized.  Note  that  the  location  of  the 
maximum  exposure  follows  a  fairly  regular  path.  The  highest 
counts  are  observed  at  the  outlet  of  the  northern  part  of 
the  valley  or  at  the  entrance  to  the  Garmisch  basin,  re¬ 
spectively. 

The  meteorological  conditions  may  be  seen  from  Figs.  A1  - 
A3  in  the  appendix  A.  Fig.  A1  shows  the  wind  speed  pro¬ 
files  (obtained  by  pibal  tracking  at  Farchant) ,  Fig.  A2 
shows  both  the  trajectories  and^  the  time-height  curves 
of  the  balloon  flights  (the  ascents  B,  C,  D  were  made 
during  the  release  from  11.30  to  12.30  CET,  see  Table  4). 

It  is  to  be  seen  that  the  up-valley  wind  was  well  de¬ 
veloped  during  the  time  of  the  experiment  and  that  a  mean 
wind  speed  of  U  s  5  m/s  is  a  reasonable  approximation  to 
the  layer  between  the  valley  floor  and  the  source.  The 
temperature  structure  (Fig.  A3)  is  characterized  by  an 
elevated  temperature  inversion  with  its  base  between  300 
and  400  m  above  the  valley  floor  (ascent  B,  D) ,  i.e.,  the 
inversion  barrier  is  present  immediately  above  the  source. 
With  overcast  conditions  a  D-stability  class  was  indica¬ 
ted. 
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The  results  are  summarized  in  Fig.  9,  which  shows  the 
ground-level  normalized  concentration  SU/Q  as  a  function 
of  distance.  The  theoretical  curves  are  centerline  values. 
The  basic  model  input  parameters  are  listed  in  the  legend 
of  Fig.  9.  Case  a  is  the  "normal"  Gaussian  plume  model 
with  total  reflection  at  the  ground;  case  b  assumes  total 
reflection  at  the  inversion  lid.  Case  c  describes  the  re¬ 
flection  from  both  the  ground  and  the  inversion,  i.e. 
the  trapping  of  the  tracers  within  this  layer;  case  d 
considers  also  the  reflection  from  the  walls  (diffusion 
in  the  channel  of  width  B  and  height  h) .  Since  the  val¬ 
ley  width  B  =  2000  m  is  relatively  large  and  the  plume 
width  keeps  relatively  small  under  D-stability  con¬ 
ditions,  there  is  practically  no  difference  between  the 
curves  d  and  c. 

The  best  approximation  to  the  measured  concentration 
pattern  is  curve  b,  assuming  perfect  reflection  from 
the  inversion,  but  none  from  the  ground.  This  does  not 
mean,  however,  that  the  diffusion  really  took  place  in 
this  way.  It  was  pointed  out  by  Fosberg  et  al.  [4]  that 
the  dispersion  may  not  only  be  restricted  to  small  scale 
turbulence  but  is  also  due  to  organized  divergence  fields 
occurring  within  the  mesoscale  valley  wind  circulation. 

The  authors  propose  a  'divergence  correction'  to  be 
applied  to  the  Gaussian  plume  model;  for  realistic  esti¬ 
mates  of  the  'toposcale'  divergences  it  was  shown  that 
this  term  would  reduce  the  concentration  maximum  by  a 
factor  of  more  than  2.  This  is  the  order  of  magnitude 
observed  in  our  field  experiment.  Another  difficulty 
stems  from  the  fact,  that  the  tracer  material  was  not 
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OISTANCE,  Km 

Fig.  9: 

Wind-speed-normalized  relative  concentration  SU/Q  as  a  func¬ 
tion  of  distance  X.  Theoretical  curves  (axial  values)  within 
the  frame  of  the  extended  Gaussian  plume  model.  Model  input 
parameters:  h  =  375  m,  B  =  2000  m,  H  =  300  m,  ZQ  =  0.40m, 
a  =  b  =  1000  m,  D  stability. 
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released  from  a  stack  but  from  the  top  of  a  hill  (acting 
like  an  obstacle  to  the  valley  wind).  Thus,  near  the 
source  wake  effects  -  with  particle  loss  -  may  occur, 
the  relatively  high  concentration  at  x  =  1,8  km  is  pre¬ 
sumably  caused  by  this  effect.  Due  to  the  distortion  of 
the  windfield  by  the  'release  mountain'  the  effective 
height  H  is  not  well  defined  too;  therefore ,  the  assump¬ 
tion  H  =  300  m  (i.e.  the  equalizing  of  the  effective 
source  height  with  the  height  of  the  hill)  is  certainly 
only  a  first  (crude)  approximation. 

The  main  facts  are,  that  the  observed  concentration  field 
shows  an  increase  with  distance  as  to  be  expected  under 
D-stability  conditions  for  the  considered  distance  range, 
and  that  the  maximum  concentration  is  within  a  factor  of 
2  in  agreement  with  the  simple  Gaussian  plume  model. 

5.5.2.  Unstable  Lapse-rate  Conditions 


Most  experiments  were  conducted  during  unstable  lapse- 
rate  conditions  with  no  inversion  structures  being  pre¬ 
sent  in  the  lower  layers.  Fig.  10  shows  the  normalized 
centerline  concentration  as  a  function  of  downwind  dis¬ 
tance  for  the  cases  with  strong  insolation  and  with  well 
developed  up-valley  winds  (see  also  Table  4) .  The  mea¬ 
sured  values  are  in  good  agreement  with  the  Gaussian 
absorption  model  (curves  for  the  B-  and  C-stabilities) . 
Compared  to  the  Gaussian  plume  model  with  total  reflec¬ 
tion  from  the  ground  the  agreement  is  again  within  a 
factor  of  2. 


3 


DISTANCE,  Km 

Fig.  10:  Centerline  ground-level  normalized  concentration 
versus  downwind  distance  for  experiments  with 
strorg  insolation  and  well  developed  up-valley 
wind.  Curves  (for  B-  and  C-stabilities)  accor¬ 
ding  to  the  Gaussian  absorption  formula  (Eqn.7b) 
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12:  Measured  ground-level  normalized  concentration 

as  a  function  of  downwind  distance  for  different 
stability  conditions. 


Fig,  1 1  shows  the  results  of  three  experiments  whose 
stability  conditions  were  ranked  between  C  and  D.  For 
comparison  the  Gaussian  absorption  model  curves  are  also 
shown  for  both  stability  classes.  There  is  no  obvious 
agreement  with  the  theory,  but  the  tendency  of  the  func¬ 
tional  behaviour  is  evident:  in  2  cases  there  is  an  in¬ 
crease  of  SU/Q  with  distance  (May  and  August  experiments) , 
and  in  the  remaining  case  there  is  a  definite  decrease 
with  distance  (as  to  be  expected  for  C-stability  condi¬ 
tions)  . 

The  range  of  the  observed  normalized  centerline  concen¬ 
trations  for  different  stability  conditions  is  to  be 
seen  from  Fig.  12.  As  already  mentioned  the  highest  SU/Q 
values  were  observed  under  D-stability  conditions  with 
an  inversion  layer  aloft.  At  the  distance  of  9  km  (i.e. 

at  the  entrance  to  the  Garmisch  basin)  a  value  of 
—6  —2 

SU/Q  =10  m  was  measured,  exceeding  those  of  the 
B-stability  cases  by  an  order  of  magnitude. 

5.5.3.  Special  Aspects 


5 . 5 . 3 . 1 .  Source  Above  a  Lifted  Ground  Based  inversion 


A  very  interesting  experiment  (No.  10)  took  place  on 
16  December  1975.  In  this  case  the  aerosol  source  was 
located  above  a  persistent  lifted  ground  based  inversion 
with  a  temperature  increase  of  about  10°C  between  100 
and  300  m  (Fig.  13) . 
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Fig.  13:  Temperatur  profiles  on  16  December  1975 
(Experiment  No.  10) . 


In  the  first  100  m  a  weak  wind  with  northern  (up-valley) 
directions  was  prevailing;  above  this  shallow  cold  layer 
a  warm  foehn  current  was  blowing  into  the  opposite  di¬ 
rection  (see  Fig.  A4  -  ascents  B,  E,  F  -  in  the  appen¬ 
dix  A) .  It  is  remarkable  that  even  then  tracers  could 
be  counted  in  the  up-valley  area  (Fig.  14) .  Obviously, 
a  small  fraction  of  the  particles  was  able  to  penetrate 
the  strong  inversion  layer  and  could  reach  the  up-valley 
flow  in  the  bottom  cold  layer. 


14:  Particle  counts  and  plume  centerline 
Experiment  No.  10  (16  December  1975) 
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5. 5. 3. 2.  Path  of  the  Plume  Centerline 

It  was  shown  in  section  5.5.2.  that  similar  meteorologi¬ 
cal  conditions  (well  developed  up-valley  wind,  strong 
insolation,  B-stability)  lead  to  comparable  normalized 
ground- level  centerline  concentrations  along  the  valley. 
This  statement  is,  however,  not  valid  for  the  whole  cross- 
section  of  the  valley,  since  the  path  of  the  centerline 
may  very  well  differ  from  one  experiment  to  the  other. 

This  is  clearly  to  be  seen  from  a  comparison  of  the  par¬ 
ticle  counts  on  28  June  1976  (Fig.  15)  and  on  7  July  1976 
(Fig.  16) .  An  inspection  of  the  counts  at  the  Farchant 
cross-section  indicates  that  in  both  cases  the  plume 
centerline  was  not  constrained  to  the  center  of  the  val¬ 
ley:  in  case  of  the  June  experiment  it  was  found  at  the 
Kramer  valley  side  and  in  the  other  case  at  the  Wank 
valley  side.  Note  that  at  the  same  time  the  particle 
counts  were  almost  identical. 

The  different  courses  of  the  particle  cloud  could  also 
be  identified  in  the  Garmisch  basin,  the  respective 
"homologous"  distribution  of  the  counts  is  evident  from 
Figs.  15  and  16. 

Note  that  in  both  cases  the  up-valley  wind  was  similarly 
well  developed  and  that  no  obvious  differences  in  the 
flow  behaviour  could  be  stated.  Nevertheless,  the  dif¬ 
ferent  paths  of  the  particle  cloud  hint  at  amesoscale 
meandering  of  the  mean  flow,  but  for  the  manifestation 
of  the  meandering  more  elaborate  measuring  techniques  are 
required . 


Particle  counts  and  plume  centerline  for 
Experiment  No.  13  (28  June  1976) 


lij  Particle  counts  and  plume  centerline  for 
Experiment  No.  14  (7  July  1976) 
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5 . 5 . 3 . 3 .  Particle  Concentrations  at  the  Mountain  Stations 


Another  important  aspect  follows  from  a  comparison  of 
the  particle  counts  at  the  mountain  stations  for  different 
synoptic  (gradient)  flow  conditions. 

On  28  June  1976  (Fig.  15)  relatively  high  particle  counts 

were  observed  at  the  Kreuzeck  area,  these  were  even 

higher  than  those  at  the  Garmisch  basin  (the  ratio  of 

the  counts  for  the  plume  centerline  amounts  to 

R  z  300/100  =3).  In  this  case  the  gradient  wind  was 

blowing  from  the  eastern  sector  and,  therefore,  was  not 

opposed  to  the  up-valley  (NE)  wind  channelled  in  the 

valley  (see  ascent  E  in  Fig.  A6) .  It  can  be  seen  from 

the  trajectories  of  the  pilot  balloons  near  the  entrance 

to  the  Garmisch  basin  (Fig.  A5)  that  the  flow  in  the 

# 

considered  height  range  was  heading  directly  towards 
the  Kreuzeck  mountain  chain,  and  this  fact  is  thought 
to  be  the  main  reason  for  the  observed  high  particle 
concentrations . 

On  7  July  1976,  however,  practically  no  particles  could 
reach  the  mountain  samplers  (Fig.  16).  In  this  case  the 
gradient  wind  was  blowing  from  the  western  sector  and, 
therefore,  had  a  component  opposite  to  the  up-valley 
wind  direction  (see  ascent  B  in  Fig.  A8) .  The  reversal 
to  westerly  directions  took  place  in  the  height  range 
of  the  mountain  stations  (see  also  the  pilot  balloon 
trajectories  in  Fig.  A7) ,  and  this  change  of  direction 
-  connected  with  low  wind  speeds  -  seems  to  be  respon¬ 
sible  for  the  absence  of  the  particles. 
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The  concentration  pattern  at  the  Kreuzeck  area  obviously 
depends  to  a  high  degree  on  the  inlet  of  the  valley  wind 
into  the  Garmisch  basin,  and  this  inlet  may  be  influen¬ 
ced  at  upper  levels  by  the  superimposed  synoptic  flow. 
This  comparison  confirms  also  the  well-known  fact  that 
the  valley  wind  system  at  low  levels  is  often  complete¬ 
ly  decoupled  from  the  upper  winds. 


6.  CONCLUSIONS  AND  FINAL  REMARKS 


The  analysis  of  the  tracer  experiments  described  here 
illustrates  the  complexity  of  pollutant  transport  in  a 
mountain  valley  but  also  demonstrates  the  difficulties 
in  separating  the  different  effective  dispersal  processes. 
The  transport  is  governed  to  a  high  degree  by  valley  and 
slope  wind  circulations  with  specific  flow  structures, 
by  enhanced  turbulence  over  rugged  terrain  and  by  along- 
valley  channelling  effects.  Since  the  interaction  of 
these  processes  is  nonlinear,  the  measured  concentration 
patterns  are  difficult  to  relate  to  model  conceptions. 

Most  experiments  were  conducted  with  well-developed  up- 
valley  wind  and  unstable  lapse-rate  conditions  (B-  and 
C-stabilities) .  It  is  found  that  the  normalized  axial 
concentration  of  the  tracer  along  the  valley  (SU/Q)  is 
in  good  agreement  with  the  Gaussian  formula  without  re¬ 
flection  from  the  surface  (Eqn.  7b)  and  with  a  sigma 
scheme  applying  to  flat  terrain  conditions.  On  the  other 
hand,  compared  to  the  classical  Gaussian  plume  model 
(where  the  ground  is  assumed  to  be  a  perfect  reflector 


* 
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of  the  tracers,  Eqn.  7a)  the  agreement  is  within  a  factor 
of  2.  Dispersion  by  bulk  fluid  motion  may  be  partly  re¬ 
sponsible  for  this  difference.  Indeed,  Fosberg  et  al. 

[4]  have  shown  that  "mass  divergence  influences  on  plume 
dispersion  may  modify  classic  Gaussian  calculations  by  as 
much  as  a  factor  of  two".  It  must  be  mentioned,  however, 
that  in  our  experiments  the  tracers  were  released  from 
the  top  of  an  isolated  hill.  Due  to  distortions  of  the 
windfield  by  the  hill  the  effective  source  height  is  not 
well  defined.  Also,  wake  effects  -  with  particle  losses  - 
are  to  be  expected  near  the  source.  Both  effects  could 
cause  reductions  of  ground-level  concentration  at  greater 
downwind  distances. 

The  highest  centerline  concentrations  of  SU/Q  (which 
exceed  those  of  the  B-stability  cases  by  an  order  of 
magnitude)  were  observed  under  D-stability  conditions 
with  an  inversion  layer  immediately  above  the  source. 

In  this  case,  a  distinct  increase  of  the  concentration 
with  downwind  distance  was  found;  this  indicates  that, 
as  expected,  the  inversion  restricted  the  vertical  dis¬ 
persion  of  the  tracer.  The  measured  concentration  pattern 
is  in  good  agreement  with  a  Gaussian  plume  model  assuming 
perfect  reflection  from  the  inversion,  but  none  from  the 
ground  (Eqn.  6). 

The  complexity  of  the  aerosol  transport  becomes  most 
evident  in  those  cases  when  atmospheric  conditions  were 
such  that  none  of  the  Gaussian  plume  models  referred  to 
above  could  be  applied.  For  example,  the  tracers  were 
able  to  partly  penetrate  a  very  strong  inversion  layer 


and  vertical  wind  shear  enhanced  this  penetration.  Se¬ 
condly,  the  location  of  the  centerline  differed  from  one 
experiment  to  another;  mesoscale  meandering  of  the  basic 
flow  is  thought  to  be  responsible  for  this.  Finally,  in 
one  case,  at  least,  the  concentrations  at  the  mountain 
stations  were  considerably  higher  than  those  in  the  valley. 
This  occurrence  appears  to  be  dependent  on  the  manner  in 
which  the  valley  wind  enters  the  Garmisch  basin  at  upper 
levels  and  this,  in  turn,  depends  upon  the  synoptic  flow 
pattern . 

More  tracer  experiments  in  the  Loisach  river  valley  area 
are  planned  for  the  near  future.  With  the  tracer  source 
at  ground  level,  the  main  interest  will  be  in  the  aerosol 
transport  in  nighttime  diabatic  mountain  winds.  The 
"straight  channel"  model  reviewed  in  the  first  part  of 
this  report  will  provide  a  theorotical  basis  for  inter¬ 
pretation  of  the  new  experimental  data. 
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APPENDIX  A 


Figures  A1  -  A8 
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Fig .  A1 ;  Wind  speed  profiles  on  8  March  1976  (Experiment  No.  12) 
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Temperature  profiles  on  8  March  1976 
(Experiment  No.  12) 
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APPENDIX  B 


Table  B1 :  Diffusive  function  F * (a* ,  1*,  £*) 

to  Eqn .  (18). 


S/SLim  =  F*(oy*'  y*»  a*)  •  F*(az*,  z*,  H*) . 


•  V .  <  v . 


according 
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ZETA*  •  O.fl  (  u*  UPPER  SCALE  ). 
L*  t  0.0  0.1  0.2  0.3  0.4 


1.0  0.9  0.8  0.7  0.6 


SIGma*  t 

0.01 

79.79 

0.00 

0.00 

0.00 

0.00 

0.02 

39.89 

0.00 

0.00 

0.00 

0.00 

0.03 

26.60 

0.10 

0.00 

0.00 

Q.QQ 

0.04. 

19.95 

0.88 

0.00 

0.00 

Q  .00 

0.05 

15.96 

•2.16 

0.01 

0.00 

0.00 

Q .  06 

13.30 

3.32 

0.05 

0.00 

0.00 

0.07 

11.40 

4.11 

0. 19 

0.00 

0 .00 

0.08 

9.97 

4.57 

0.44 

0.01 

0.00 

0.09 

8.87 

4.78 

0.73 

0.03 

Q.OQ 

0.10 

7.98 

4.84 

1.08 

0.09 

O.OQ 

0.11 

7.25 

4.30 

1.39 

0.18 

0.01 

0.12 

6.65 

4. 70 

1.66 

0.29 

0.03 

0.13 

6.14 

4.57 
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). 
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